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In the late stages of nuclear burning for massive stars (M > 8 Mq), the 
production of neutrino-antineutrino pairs through various processes becomes the 
dominant stellar cooling mechanism. As the star evolves, the energy of these 
neutrinos increases and in the days preceding the supernova a signihcant fraction 
of emitted electron anti-neutrinos exceeds the energy threshold for inverse beta 
decay on free hydrogen. This is the golden channel for liquid scintillator detectors 
because the coincidence signature allows for signihcant reductions in background 
signals. We hnd that the kiloton-scale liquid scintillator detector KamLAND 
can detect these pre-supernova neutrinos from a star with a mass of 25 Mq at a 
distance less than 690 pc with 3a signihcance before the supernova. This limit 
is dependent on the neutrino mass ordering and background levels. KamLAND 
takes data continuously and can provide a supernova alert to the community. 

Subject headings: neutrinos, supernovae: general 


1. Introduction 


The hr st extr asolar neutrinos wer e detected fror n SN1987A by the Ka. miokande-11 (IHirata et al 


19871. Il988 h IMB ( Bionta et al.11198711 . and Baksan ( Alekseev et al.lll987h experiments. This 
dataset h as provided rn any insights into the properties of neutrinos and the physics of su¬ 
pernovae fjVissanil 120151) . SN1987A was located in the Large Magellanic Cloud at a distance 
of ~ 50 kpc. A core-collapse supernova in the Milkyway proper would provide a larger hux 
of neutrinos. This combined with the large suite of r unning neutrino experiments makes the 
next Galactic supernova a greatly anticipated event (jScholbergi 120121) . 


In a Type II supernova, a huge burst of neutrinos is released, carrying away ~ 10^^ ergs of 
energy in 10 s. Leading up to this cataclysmic event, neutrinos have already been playing an 
important role in the cooling of the evolving giant star. Starting in the carbon burning phase, 
the dominant mechanism for cooling these massive (M > 8 Mq) stars is the los s of energy 


due to n/ n pairs created by t herm al processes. From applying the discussion in lltoh et ah 


fll996l) to IWooslev fc Hegerl (120151) . the dominant process in most of M > 10 Mq stars 
is the pair process, e’''e“ —)■ vv. For other stars with smaller masses, the plasmon decay 
becomes more important, j ^ uu. Secondary contributions come from the photo process, 
76 “ —)■ e~uu and bremsstrahlung, e~{Ze) —)■ {Ze)e~uu. These thermal processes are often 
used to set limits on non-stan dard neutrino inte ractions since such processes would change 
the evolution of these objects (iHeger et al.l 120091) . The most stringent limit s on the neutrino 
magnetic moment come from this type of analysis (lArceo-Dfaz et al.l 120151) . 
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Since these thermal neutrinos precede the supernova, they can also be called pre¬ 
supernova neutrinos (preSN). Figure [T] shows the overall time evolution of the Vf. luminosity 



though the preSN luminosity is several orders of magnitude smaller than the SN luminosity, 
the detection of preSN is desired since the preSN encode information about the late stages 
of stellar evolution for high mass stars and could act as a supernova alert, a more detailed 
discussion is found in section [5l 


SN extend to a few tens of MeV. In comparison, the average energy of preSN is low, 
typically E < 2 MeV. In this energy range, there are three reactions that can be used to 
detect these neutrinos in realtime: coherent neutrino scattering, neutrino-electron scattering 
and inverse beta decay (IBD), z/g + p ^ -|- n. IBD has one of the highest cross sections 

for neutrino detection. It also has relatively low backgrounds due to the easily identihable 
delayed coincidence signal created by the prompt positron annihilation followed by the de¬ 
layed neutron capture. Depending on the detector material, coherent neutrino scattering 
may have a higher cross section than IBD, but the signal has never been observed due to 
the very low reconstructed energy of the recoiling nucleus. The detection of preSN through 
neutrino-electron scattering is possible. However, its cross section is lower than IBD, which 
reduces the total number of detected events, and the background rate is high since there is 
no coincidence signal. Thus, IBD is the most promising channel for preSN detection. 


The energy threshold for IBD is 1.8 MeV. A few days before the supernova, a signihcant 
fraction of he exceeds the IBD threshold and it becomes possible to detect the preSN with 
IBD. IBD is the main supernova channel for both liquid-scintillator detectors and water- 
Cherenkov detectors like Super-Kamiokande. Water-Cherenkov de tectors have relatively 
high energy thresholds, such as Ef, = 4.5 MeV flRenshaw et al.ll2014l) . This limits both the 
number of IBD prompt events and the efficiency for detecting the delayed neutron capture. 
In comparison, monolithic liquid-scintillator detectors have energy thresholds below 1 MeV 
and are therefore able to sample a larger fraction of the preSN prompt energy spectrum and 
effectively detect the neutron capture. Thus, liquid-scintillator detectors have an advantage 
for the detection of preSN, even if they are smaller than typical water-Cherenkov detectors. 


There are two operating mo nolithic liquid-scintilla tor detectors with low-en ergy thresh- 
olds, KamLAND and Borexino ( Cadonati et ah 2002 1. The SNO-I- detector f Chen 2008) 
is e xpected t o come online soon and construction has started on the 20 kton JUNO detec¬ 
tor (iLil 1201411. In add ition, there are several proposals for multi-kt on experiments su ch as 


RENO-50 f Kimll2014ll . HANOHANO ( Learned et al.ll2008ll . LENA f Wurm et al.ll2012[l . and 
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ASDC flAlonso et al.ll2014l) . All of these detectors would be sensitive to this preSN IB P sig¬ 


nal. A large Gd-doped water-Cherenkov detector such as Gd-doped Super-Kamiokande flBeacom fc Vagins 
20041) would have increased sensitivity due to the higher neutron capture detection efficiency 
but the higher energy threshold continues to limit the sensitivity. The Baksan and LVD 
scintillator detectors are similarly limited in their sensitivity to preSN d ue to their relatively 
high energy thresholds (iNovoseltseva et ahl 1201 ll: lAgafonova et al.l 120151) . 


In the previous studies fjOdrzvwolek et al.ll2004 lOdrzvwolek fc Hegeiil2010l: iKato et ah 


20151 ). the expected number of IBP events in several detectors is evaluated without a detailed 
detector response model. We focus on KamLANP since it is currently the largest monolithic 
liquid-scintillator detector. In this article, we quantify KamLANP’s sensitivity to the preSN 
using the actual background rates and a realistic detector response model. We discuss 
the development of a supernova alert b ased on preSN. Betelgeuse is a well-known possible 
supernova progenitor (IPolan et al.ll2014i) and we evaluate the performance of the preSN alert 
based on this astrophysical object. 


2. PreSN signal 


The hrst calculation of the n umber of detected preSN is foun d in lOdrzvwolek et al. 


(120041) and updates can be found in lOdrzvwolek fc Hegerl fl2010l) and iKato et al.l (120151) . We 
use the preSN spectra 0 m(^, d) as a function of time and energy from the Odrzywolek’s 
results corrected for the distance of d to the pre-supernova star. We use this to calculate 
KamLANP’s sensitivity to preSN with two example stars of M = 15 Mq and M = 25 Mq. 
Figure [2] shows the time evolution of the i>e luminosity in the top panel and the averaged 
energy in the middle panel during the 48 hr before the collapse. The integrated luminosity 
over the last 48 hr preceding collapse is 1.9 x 10^° erg and 6.1 x 10®° erg, respectively, for 
the two star masses. They correspond to 1.2 x 10®° and 3.8 x 10®° Ve, respectively. The 
weighed differential luminosity by energy, Ep^dL/dEp^ ~ dL/dlogEp^, is also shown in the 
bottom of Figure [2] with the SN for reference. The average energy of the integrated Ue hux 
is 1.4 MeV and 1.2 MeV for the 15 Mq and 25 Mq models, respectively. 

In detectors, the reconstructed prompt (positron in IBD) spectrum can be written as, 
d‘^NM{t,E^^^;d) 


dtdEl^^ 


= 6n.e6s(i?r)^T / Ep^' d)R{El^^, E^^^)a{Ep^), (1) 


where Ny is the number of targe t protons in the analysis volume, a{EpJ is the IBD cross 
section (IStrumia fc Vissanil 120031 ). euve is the mean livetime-to-runtime ratio, is the 


total detection efficiency. The details of these parameters are presented in Section [3l E^^ 


IS 













































5 


the reconstructed energy and — 0.78 MeV is the expected energy of the prompt 

event for an input z/g with an energy of Ep^. The integration in Equation ([1]) is a convolution 
of the theoretical spectrum (pM with the detector response. We model the detector response 
as a Gaussian with R{E^^, E^^^) with energy resolution cte: 




\/^(Te(E) 


exp< 


(^rec _ 

2ai(E) 


( 2 ) 


Assuming preSN emitted from the 200 pc star and the perfect 1 kt detector with euve = 1, 
€s(Ep°^) = 1, and R(Ep^^, = S(Ep^^ — E^p), the event spectrum integrated over the last 

48 hr preceding collapse is shown in the left panel of Figure [3l The number of arrival z/g in 
the last 48 hr is 1.8 x 10^^ for the 15 Mq star and 6.5 x 10^^ for the 25 Mq star, respectively. 
The total number of preSN events in the detector is 44 and 95, respectively. 


The right panel of Figure [3] shows the preSN spectrum integrating over the last 48 hr and 
SN spectrum integrating over 10 sec with the vertical axis of ~ dN/d log E^'^'^. 

Our assumption of supernovae at 200 pc is 50 times closer than the usually assumed 10 kpc. 
We note that the SN from 200 pc supernovae will create ~ 10® events in the detector. The 
current KamLAND electronics will not be able to record more than the basic hit information 
for the SN. The information from the preSN will not be lost in the case of the DAQ crash. 


Neutrino oscillation in the preSN emission region redu ces the i^e flux. The flux of 
observable z/g can be expressed following iKneller et al.l (120081) as. 


'>p. =p4>t + (1 -p)0L, 


(3) 


where 0®^ and 0®^ are the original spectra of z/g and G' t . With the assurnption of 0®^ = 
0.190®^ based on 0°^ oc 0®^ and a Ux/^e ratio of 0.19 from lOdrzvwolek et al.l (1200411 . and an 
adiabatic a pproximation for p w ith sin^ 612 = 3.08 X 10“^ and sin^^ig = 2.34 x lO'^ (2.40 x 
10“^) from ICapozzi et al.l (1201411 . we have 0p^ = a0p^, where a = 0.74 (0.21) corresponding 
to the normal (inverted) neutrino mass order. The corrected spectrum for the neutrino 
oscillation is then given by multiplying the coefficient a. With the integral over the last 
48 hr preceding collapse, the event spectrum is shown with the oscillation effect and the full 
detector response in the middle and bottom panels of Figure 0] with KamLAND using the 
parameters describe in Sec. 01 


3. KamLAND detector and its background 

KamLAND is located in the Kamioka Mine in Japan’s Gifu prefecture (36.42°N, 137.31°E). 
Mt. Ikenoyama rises ~ 1 km above the detector reducing backgrounds due to cosmic rays by 
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five orders of magnitude. The KamLAND detector consists of approximately 1 kt of liquid 
scintillator, a mixture of 20% psuedocume and 80% dodecane. It is contained in a 13-m- 
diameter spherical balloon made of a 135-/im-thick transparent nylon ethylene vinyl alcohol 
copolymer (EVOH) composite him. An array of photomultiplier tubes (PMTs) is used to 
detect the scintillation light from events occurring within the balloon. This array consists of 
1,325 fast PMTs masked to 17-inch diameter to achieve the desired timing performance and 
554 older 20-inch diameter PMTs reused from the Kamiokande experiment. The PMTs are 
mounted on the inner surface of an 18 m-diameter stainless steel sphere. Non-scintillating 
mineral oil hlls the space between the balloon and the inner surface of the sphere. Its den¬ 
sity is matched to the liquid scintillator to support the balloon and also acts as a passive 
shielding against external backgrounds from the sphere, PMTs, and surrounding rocks. This 
inner detector is further shielded by a 3.2kton water-Cherenkov veto detector. In 2011, a 
3.08 m-diameters inner balloon containing 13 tons of Xe-loaded liquid scintillator (Xe-LS) 
was installed in the center of the main balloon as a part of the Karn LAND Zero-Neutrino 
Double-Beta Decay (KamLAND-Zen) experiment (IGando et al.l 120121) . 


The position and energy of an event within the balloon can be reconstructed using the 
timing and charge distribution obtained from the PMT array. The reconstruction is cali¬ 
brated by a number of radioactive sources: ^°C o, ®®Ge, ®^Zn, ^^^Am^Be, ^^^Gs, and 

210pq13q flBerger et al.ll2009l: iBanks et al.ll2015l) . From these calibrations and naturally oc¬ 
curring radioactive sources, the energy resolution (crE(F'))is determined to be 6.4 %/\/E (MeV) 
and the position reconstruction 12 cm/y^F^^MeV). 

Gandidate Ve events corresponding to the prompt positron annihilation and delayed 
neutron capture of the IBD interaction are selected with a series of cuts on the energy, 
position, time and space coincidence of the two events. The two events must occur with 
0.5 < AT(/rs)<1000 and within Ai? <2.0 m, where AT and AT are time and spatial 
differences. The reconstructed position of both events must be within a spherical hducial 
volume Tp, Td < 6 m, which determines the hducial number of target proton. At = 5.98 x 
10^^. The reconstructed energy of the prompt event is required to be in the energy range 
T'p®‘^(MeV) > 0.9. The delayed event has an energy characterized by the energy of the 
neutron capture gammas. Two energies are used: 1.8 < T'^®'^(MeV) < 2.6 corresponding to 
capture on H and 4.4 < T'^®'^(MeV) < 5.6 corresponding to capture on ^^G. An energy and 
position dependent likelihood variable is constructed to differentiate Vg from backgrounds 
due to accidental coincidences, w hich become more likely at lower energies and as events 
reconstruct closer to the balloon flGando et al.N2nill ). Finally, an additional position cut 
on the delayed event is applied to eliminate backgrounds due to the KamLAND-Zen inner 
balloon and support structure. The cut eliminates a central sphere and cylinder: < 2.5 m 

and \/x\ + y\ < 2.5 m for Zd > 0 m where (a;d,2/d,^d) is the reconstructed position of the 
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delayed capture event. 


This series of cuts matches the standard KamLAND analysis flGando et ahl 120131) . The 
minimum prompt energy is chosen to guarantee 100 % detection efficiency for z/g with a 
energy of 1.8 MeV. The total efficiency of these cuts, es(i7p®^), is energy dependent due to 
the Likelihood selection as shown in Figure |4] (Top). The efficiency loss is dominated by the 
inner balloon cut. Without this cut, the efficiency is higher, ~0.9 depending on the energy 
of the prompt event. Further cuts are used to remove backgrounds due to high energy muon 
events. The effect of these cuts is to reduce the effective livetime to enve = 0.903 even though 
KamLAND takes data continuously. 


After these cuts are applied, the measured event spectrum in the preSN region (0.9 < 
F'p®'^(MeV) < 3.5) is mainly from reactor W and geological W produced in the Earth’s interior. 
These are the backgrounds to the preSN signal. Since the Great East Japan Earthquake, the 
reactors in Japan have been off. This is the low-reactor phase, and the reduced backgrounds 
increase our preSN sensitivity. The measured background spectra are shown in the middle 
and bottom panels in Figure 0] for the low-reactor phase and the high-reactor phase. In this 
hgure, the background spectra are normalized to the 48 hr window. 


4. Sensitivity 

The middle and bottom panels in Figure |4] show measured background and expected 
preSN spectra, integrated over the 48 hr window immediately before the collapse. Differ¬ 
ent integration lengths were studied and this length was chosen to maximize the signal-to- 
background ratio. To study the sensitivity, we use the analysis range of 0.9 < E'p®'’(MeV) < 
3.5 to maximize the preSN signal to background ratio while retaining >90% of the preSN 
signal. The background rate is i?iow = 0.071 event/day in the low-reactor phase. If the reac¬ 
tors return to normal operations, the background rate rises to i?high = 0.355 events/day. The 
average efficiency is (eg) = 0.64 in this range. For the number of events N^{d) = aNM{d), 
where Z indicates the neutrino mass ordering (Z = N/I for the normal/inverted order), 
we integrate Eq. ([1]) from 0.9 MeV to 3.5 MeV and over the 48 hr before the collapse. 
They are now iV/^jy^^(200 pc) = 12.0, ^^"^^^(200 pc) = 25.7 and iV| 5 jy ^^(200 pc) = 3.38, 
^ 25 ^ 0 (200 pc) = 7.28. The Nj^i^d) are shown as a function of d in Figure [5l 

Using Nj^i^d) with background 2 daysxiJj (i=low or high), the corresponding detection 
significance, expressed in sigma, is estimated. The results for a detection significance of 
Icr, 3a, and 5a assuming i?iow are plotted in Figure |5l If a 3(T significance is required for 
preSN detection, KamLAND is sensitive to preSN from a 25 Mq star at 690 pc assuming 





Blow and the normal neutrino mass ordering. In the worst case, KamLAND is still sensitive 
to a 15 Mq star at 250 pc. In this distance range, there are several red-supergiants which 
could lead to supernovae: Antares (150 pc), Betelgeuse (200 pc). Epsilon Pegasi (210 pc). 
Pi Puppis (250 pc). Sigma Canis Majoris (340 pc), NS Puppis (520 pc), CE Tauri (550 pc) 
and 3 Ceti (640 pc). 


Be telgeuse has been stud i ed extensively as a nea rby pre-supernova star, see lTownes et ah 


(I 2 OO 9 I) ; iHaubois et al.l (120091) ; lOhnaka et al.l (120091) , therefore we use it to determine Kam- 
LAND’s sensitivity as a function of time before collapse. Bete lgeuse’s measured mass 
M = 17-25 Mq and distance d = 197 ±45 pc are highly correlated (iHarper et al.ll2008l) . We 
studied the two extreme cases: (M, d) = (15 Mq, 150 pc) and (25 Mq, 250 pc). The ex¬ 
pected time evolution of signihcance with the 48 hr integration window is shown in Figure El 
assuming the low-reactor background. 


If Betelgeuse has (M, d) =(15 Mq, 150 pc), KamLAND will easily detect its preSN. 
A 3a detection of preSN would be 89.6 (7.41) hr before collapse for the normal (inverted) 
mass ordering. If Betelgeuse has (25 Mq, 250 pc), the increased distance reduces the preSN 
flux and the number of hours before collapse which KamLAND could detect preSN. If the 
reactors in Japan are restarted, the number of hours is also reduced because of the larger 
backgrounds. Tabled] summarizes the results and shows for all of these cases that KamLAND 
can still detect preSN and has the ability to send a supernova alarm before collapse. 


5. Discussion 


We quantitatively evaluated KamLAND’s sensitivity to preSN. Using Betelgeuse as a 
likely progenitor, we hnd that a 3a detection of preSN at 2-90 hr before the collapse is 
possible. This a great improvement over the warnin g time provided by the current Super- 


Nova Early Warning System (SNEWS) described in lAntonioli et al.l (120041) . The SNEWS 


alarm is not sufficient to bring gravitatio nal wave and neutrino detectors back online from 
commissioning or calib ration campaigns JPoolev 2015 ) or redirect telescope time for shock 
breakout observations ( Adams et ah 20131) . In comparison, the 2 hour to days provided by 
KamLAND’s detection of preSN would facilitate these measurements. 


This is the motivation for the development of the preSN alert system. It provides two 
levels of alarms. A low-level alarm is produced using the semi-realtime preSN detection 
signihcance. This is calculated using a 48 hr integration window and the background level 
averaged over the past three months. A new window is opened every 15 min and has a 
latency of 25 min due to KamLAND’s online data processing time. This system requires 
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users to sign up to receive the current significance of detection, see the KamLAND web site0. 
A high-level alarm will report any 5a detections to the The Gamma-ray Coordinates Net¬ 
work (GCN) and/or the Astronomer’s Telegram (ATel). This alarm is only sent after that 
the collaboration rejects other possibilities such as DAQ problems, mine activity, and radon 
contamination. Unfortunately, the IBD reaction does not provide directional information so 
a definitive localization requires SNEWS alarms, direction detection in Super-Kamiokande, 
coherent network analysis of gravitational waves (IHavama et ahl 120151 ) and/or electromag¬ 
netic observations. However, the limited number of known local progenitors limits the total 
number of targets that could create a preSN signal in KamLAND and therefore some local¬ 
ization is still possible. 

In addition to the alarm, preSN could teach us about neutrino and supernova properties. 
The detection of preSN from Betelgeuse’s supernova and an improvement in the measurement 
of its mass and distance would allow the determination of the normal neutrino mass ordering 
at 3.6(7 (2.5(7) for the 25 Mq (15 Mq) star. The dete ction of preSN co uld also distinguish 
between a supernova with an ONe core or an Fe core flKato et al.l 120151) . 


6. Summary 


We evaluated the KamLAND’s sensitivity to preSN, neutrino-antineutrino pairs from 
massive stars that have entered the late burning phases. Detection of preSN could provide 
an early warning for the imminent occurrence of a supernova and provide an opportunity to 
study stellar evolution models in the last stages prior to collapse. This study quantitatively 
confirms that the 3cr detection of preSN is possible for stars with distances up to 690 pc 
under optimal conditions. The number of pre-supernova stars at this distance is limited but 
includes several promising candidates such as Betelgeuse. KamLAND presently provides the 
community with a preSN alarm based on the semi-realtime significance calculation and will 
send a report to GGN/AGTel for any 5a detection that have been verified by the KamLAND 
collaboration. 


There are several oth er detectors under c onstruction or pr oposed: S NO-I- flGhen 


RENO- 50, HANOHANO ( 


^earned et aLll2008l) . JUNO (lLill2014l) . L ENA (Wurm et ah 


ASDG flAlonso et al.! l2014l) . and Gd-doped Super-Kamiokande (IBeacom fc Vagins 


2008), 


2012 ), 


2004h . 


These detectors are expected to have similar or higher sensitivity to preSN and a coinci¬ 
dent preSN search would significantly reduce false signals. A future combined alarm system 
could increase the detection range to a few kpc. This extended range would include several 


^http://www.awa.tohoku.ac.jp/kamland/SNmonitor/regist/index.html 
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other pre-supernova stars including Eta-Carinae. 

The preSN are an exciting tool for both the study of stellar evolution and the supernova 
observation. The current models of preSN production are fairly simple. KamLAND’s sensi¬ 
tivity to both the flux and energy of the preSN could be used to extract more information 
on the late stages of evolution if more detailed predictions become available. 
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Fig. 1.— Time evolution of the Ve luminosity of preSN just befo re col¬ 

lapse (jOdrzvwolek fc Hegerl 120101) and of SN after collapse (INakazato et al.l 120131 ). Note, 
the time scale of the horizontal axis which is linear after the collapse but logarithmic before 
collapse. 
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Fig. 2.— Details of preSN with 15 Mq and 25 Mq star models. (Top) time evolution of 
Vf. luminosity, (Middle) time evolution of the averaged energy, and (Bottom) differential 
luminosity weighted by energy {Ep^dL/dEp^ ~ dL/dlogEpJ integrating over 48 hr preceding 
collapse with SN integrating over 10 s for reference. SN is calculated from a Fermi-Dirac 
distribution of the flux with a total luminosity of 5 x 10^^ erg and an average energy of 
15 MeV. 
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Fig. 3.— PreSN z/g event spectrum on Earth integrating over the last 48 hr preceding the 
collapse, with assuming a distance of 200 pc and the perfect 1 kt detector. Both panels 
are basically same with the linear-scale horizontal axis (left) and the log-scale horizontal 
axis (right). The vertical axis in the right panel is weighted by to be ~ diV/dlogEp®®. 
The total number of preSN events in the detector is 44 and 95 for the 15 Mq star and the 
25 Mq at 200 pc, respectively. In the right panel, for comparison, the weighted supernova 
z/g event spectrum integrating over 10 s are also shown. Expected number of SN events is 
about 8 X 10^. The effect of neutrino oscillation is not considered. 
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Fig. 4.— (Top) The energy dependent detection efficiency es{E^^) used in the analysis. 
(Middle) The integrated preSN energy spectrum and measured background in the last 48 hr 
before collapse. The backgrounds are dominated by reactor and geological Fg. Background 
levels for the low-reactor phase (i?iow) and high-reactor phase (i?high) are shown. (Bottom) 
Same as the middle panel but shown on a log scale instead, which clearly show the background 
spectra. Our analysis range is 0.9 - 3.5 MeV in the prompt energy 
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Fig. 5.— Expected number of IBD events detected in KamLAND during the 48 hr be¬ 
fore collapse including neutrino oscillation effects as a function of distance. The flux for a 
15 Mq and 25 Mq star is shown assuming the normal and inverted neutrino mass ordering. 
Horizontal dotted lines are the signihcance of the detection (see text). 
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Fig. 6.— Time evolution of significance in the low-reactor phase. If Betelgeuse has a 
mass of 15 Mq aX d = 150 pc, the Sa detection time (shown by a dotted horizontal line) 
is 89.6 (7.41) hr with the normal (inverted) ordering before collapse. KamLAND can detect 
preSN 17.0 (4.54) hr before the collapse of Betelgeuse (25 Mq, 250 pc) at the 3a level. 
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Table 1. Expected time before the Betelgeuse supernova to reach 3cr conhdence based on 

PreSN under various parameter assumptions. 


Mass [Mq] 

Distance [pc] 

Mass ordering 

Reactor status 

Time before collapse [hr] 

15 

150 

Normal 

low 

89.6 

15 

150 

Inverted 

low 

7.41 

25 

250 

Normal 

low 

17.0 

25 

250 

Inverted 

low 

4.54 

15 

150 

Normal 

high 

46.0 

15 

150 

Inverted 

high 

3.17 

25 

250 

Normal 

high 

11.1 

25 

250 

Inverted 

high 

1.93 






